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Abstract: We investigate the effects of an all order QCD resummation of soft gluon 
emissions for Higgs boson production in association with two hard jets. We consider both 
the gluon-gluon fusion and weak boson fusion processes and show how to resum a large 
part of the leading logarithms in the jet veto scale. Our resummation improves on previous 
analyses which also aim to include the effects of multiple soft gluon radiation. In addition 
we calculate the interference between weak boson fusion and gluon-gluon fusion and find 
that it is small. 
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1. Introduction 

Searching for the Higgs boson is one of the main tasks of the LHC and its production in 
association with two low-angle (tag) jets should provide an effective way of reducing back- 
grounds. In the literature, most of the attention has been focussed upon the contribution 
from weak boson fusion (WBF), see Fig. 0(a), although the Higgs can also be produced 
through gluon fusion, as illustrated in Fig. 0(b). Detecting the Higgs boson in association 
with tag jets should provide the opportunity to gain important information on the coupling 
of the Higgs to weak bosons [1-6] and even to the top quark [7] . 

To obtain these couplings separately one would like experimentally to differentiate the 
two cases. Fortunately there are several differences which arise. In WBF, the two highest 
k± jets (which we define to be the tag jets) tend also to be the high rapidity quark jets which 
are responsible for radiating the fusing weak gauge bosons. This is to be contrasted with 
the case of gluon fusion where the highest k± jets need not be energetic low-angle quark 
jets. As a result, the tag jets in gluon fusion tend to be much less forwardly peaked and 
form a much lower dijet invariant mass than the tag jets present in WBF. The dependence 
upon the azimuthal angle between the tag jets is also sensitive to the difference between 
WBF and gluon fusion. In particular, the tag jets in the WBF channel are approximately 
independent of the azimuthal angle whereas in the gluon fusion channel they have a cos^ (j) 
dependence (for a CP even Standard Model Higgs). 

In this paper, we will focus on another important difference, namely the tendency for 
the weak (color singlet) exchange to radiate fewer particles into the region between the tag 
jets. Indeed, this difference is often exploited in order to produce a WBF rich sample of 
events; typically one insists that there be no additional jets in between the tag jets with 
the transverse momentum above some value, Qq ~ 20 GeV. The tendency of WBF to 
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Figure 1: The lowest order q + q^H + q + q diagrams in the eikonal Umit, (a) for weak boson 
fusion; (b) for gluon-gluon fusion with the Higgs being produced via a top quark loop. 

radiate less has been studied previously in the context of Higgs plus two-jet production 
in association with the emission of one extra parton [8,9], by employing parton showering 
from event generators [10], by using the truncated shower approximation (TSA) [11-16] 
and via calculations based on the soft gluon exponentiation approach [14,16]. Fixed order 
calculations [6,8,9, 17] inherently suffer from the problem that as soon as cj^j"^*+^ ^ (j2jet+H 
then an all order resummation is desirable. In this paper we shall show how to improve 
fixed order calculations by summing to all orders in the strong coupling what is likely to 
be the dominant subset of the leading logarithms in the ratio of the hard interaction scale 
and the jet veto scale Qq. Our calculations ought to improve upon the existing methods 
when it comes to accounting for multiple gluon emission in between the tag jets. Before 
we present our calculations we shall briefly review the aforementioned approaches. 

In the truncated shower approximation [11-16], fixed order QCD is used to compute 
the cross-section for the real emission of one additional parton, i.e. Hjjj. Two of the 
partons are the tag jets and the presence of the third (mini-)jet, of transverse momentum 
Pt, permits a study of the role of a jet veto in the region between the tag jets. However, 
integrating the three parton cross-section over the phase-space of the third (mini-)jet will 
lead to a divergent cross-section. In the TSA, this divergence is regulated by modifying 
the three parton cross-section (ai) to 



The scale ptsa is chosen so that the integrated cross-section is equal to the two parton 
cross-section, i.e. 



where o"o is the lowest order cross-section (either for WBF or gluon fusion). This is clearly 
a rather crude phenomenological model. It has the virtue that it contains the complete 
information of the fixed order calculation for single parton emission however it will fail 
for small enough veto scales, Qq ^ Ptsa, whereupon multiple parton emission starts to 
become important. The enhanced probability for radiation in gluon fusion events leads to 
a larger ptsa than in the case of WBF. 

The soft gluon exponentiation model [14, 16] provides for a more sophisticated treat- 
ment in the case of multiple emissions by exploiting the fact that the leading behaviour 
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comes from the emission of soft gluons which are assumed to exponentiate. The cross- 
section for emitting no gluons with transverse momenta above Qq is given by 



da = exp 




(1.3) 



where ai is the cross-section for emitting one extra gluon into the veto region. As we shall 
see, this approach docs contain an important part of the leading logarithmic soft gluon 
contribution, although it docs not include any colour mixing. 

Like the soft gluon exponentiation model, parton shower algorithms (as implemented 
in the all-purpose event generators [18, 19]) succeed in systematically improving beyond 
fixed order by including a major part of the colour diagonal soft gluon evolution. However, 
the parton showers also do not include any colour mixing contributions and the soft gluon 
evolution is usually performed in the large N approximation^. We should add that event 
generators have the virtue that they deal with hadronization, which is of course not possible 
in a purely perturbative treatment. Hadronization is, however, expected to be "gentle" in 
the sense that hard enough jets axe not much affected by it. Nevertheless, the role of 
multiple parton interactions and the soft underlying event do lead to important features 
of the hadronic final state which are inaccessible to perturbative QCD. It is common to 
chose the veto scale Qo to be sufficiently large (e.g. ^ 20 GeV) so that sensitivity to the 
underlying event is reduced [20,21]. We shall not consider the underlying event any further 
in this paper but its impact clearly needs to be accounted for when it comes to comparing 
to data. We shall also assume that the Higgs boson is sufficiently long lived that we do not 
need to consider the complexities which arise upon its decay to coloured particles. 

An alternative, experimental based, approach which should help us to understand the 
radiation pattern in Higgs production has been proposed in [16,22]. Here it is suggested to 
probe the mini-jet distribution in hard color singlet scattering events using measurements 
of the process qq ^ qq + Z via 1^-fusion. 

Soft, wide angle gluon resummations similar to those we perform here (but without 
Higgs production) have been presented in [23-28] . Generically, one is interested in restrict- 
ing radiation from a hard scattering at scale Q so that, in some specified region of the final 
state phase-space, there should be no radiation above the veto scale Qq. The goal is at 
present limited to summing the leading logarithms, i.e. terms of the form 



We shall usually find it convenient to display our results in terms of the evolution variable 



At the LHC a typical Qo is expected to be 20 GeV. Assuming a hard scale of 200 GeV 
then gives ^ ~ 0.3 with n/ = 5 and as{mz) = 0.12. If Qo can be pushed down to 5 GeV, ^ 

^A'^ is the number of colours 



:^log(Q/Qo) 



(1.4) 




(1.5) 
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could reach well above 0.5. Recall that the larger one can make ^ the more dominant one 
expects the WBF contribution to be. 

In this paper we will make an all order resummation of a part of this leading logarithmic 
series. We will, however, not include those leading logarithms which arise by virtue of the 
fact that the process we consider is "non-global" [29,30]. The problem of summing the 
leading non-global logarithms is a thorny one since the colour structure makes anything 
other than a leading N approximation impossible at the present time. Fortunately for Higgs 
production, the impact of the non-global logarithms is expected to be small provided we 
do not choose a particularly small cone radius [31, 32]. For example, in [32] the non- 
global logarithms were estimated in the two-quark case (e.g. e^e~ qq) and, for a cone 
radius R = 1 and evolution variable = 0.9, found to be a 5% correction to the primary 
emission result. For smaller ^ values the effect diminishes still further since the non-global 
logarithms enter first at order Apart from the non- global logarithms, we shall also 
neglect a leading contribution to the primary emission contribution which arises as a result 
of implementing a jet algorithm [32,33]. This too is a effect and therefore is smaller at 
lower values of ^. The results quoted in [32] indicate that these logarithms begin to play a 
role for relatively large ^ values, i.e. ^ ^ 0.5 for R = 1 (diminishing as R^ for smaller R). 
Following [32], these jet algorithm dependent corrections can be systematically computed 
and so the results we present here could be further improved upon if need be. Finally, we 
shall ignore the potential impact of the "super-leading logarithms" discussed in [34]. At 
the present time, the existence of these logarithms is uncertain. In any case, one might 
expect any such contribution to the Higgs plus jet production process to be modest since 
they enter first at order and are subleading in N. 

We will of course re-confirm that the radiation pattern is very different for WBF and 
gluon-gluon fusion. We will also show that the LHC will explore the kinematic region 
where higher order soft gluon corrections are important and fixed order calculations are 
inadequate. As for the interference (between gluon fusion and WBF contributions), this 
turns out to be surprisingly small. This is in part a result of the soft gluon evolution, 
but it is also a consequence of a tree level helicity dependent cancellation present in the 
eikonal approximation. The source of interference we consider is quite distinct from the 
crossed-channel interference discussed in [35], which vanishes in the eikonal approximation. 

The outline of this paper is as follows. First in Section || we have a quick look at the 
lowest order cross-sections to set the scene for what follows. Then, in Section |3|, we explain 
the relevant theory and show how to accommodate the extra outgoing Higgs boson. We 
also briefly discuss the recent two- loop calculations of [36]. Numerical results are then 
presented and commented upon in Section ^. Finally, we draw our conclusions in Section 

i 

2. Lowest order 

For both electroweak ZZ and WW fusion processes, shown in Fig. ||a, and the QCD gluon- 
gluon fusion process, shown in Fig. P^, we shall compute the lowest order cross-section in 
the eikonal approximation. This is a good approximation in the case where the final state 
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jets are at low angles [37]. In any case, our goal in this paper is not to provide precise 
predictions for the LHC rather it is to show how to include the leading effects of QCD 
radiation. In the eikonal approximation, only the t-channel contributions illustrated in 
Fig. 1^ are important; the vertex factor U7^n simplifies to u'y^u oc 2p^ and the parton 
helicity goes unchanged. In the electroweak vertices the factor 7^ results in an extra minus 
sign in the case of left handed fermions, i.e. ulJuI^ul — 2p^. 

Using the notation of Fig. |^ and the conventions of [38] the WBF amplitudes are then 



Mzz 

Mww = , . ... (2.1) 



(m| + k\^{w?^ + v cos^ sin^ 

2" 



(m^ + Kx. ) ("^w^ + ^2± ) ^ sin^ dw 

with V being the Higgs boson vacuum expectation value, = {GfV^)~^ ■ Here it is 
to be understood that the charged current amplitude vanishes for right-handed quarks. 
In principle the amplitude should contain CKM matrix elements but since we are not 
interested in the flavor of the outgoing quarks we neglect them. In the case of Z boson 
fusion the factors 

h = Tl,-Qism^0w (2.2) 

account for the helicity /flavour dependent factor from the Z- vertex (with T^- = for 
right-handed quarks). 

For the gluon fusion process we work in the — > 00 limit [39-41] and note that the 
approximation is good at the level of a few percent even for niH = nit [42]. It has also 
been shown to be a good approximation even in the limit where the dijet invariant mass 
is large provided k± ^ mt [42]. Up to a colour factor, the amplitude is (for quark or gluon 
external states) 

M,, = (2.3) 

Comparing the WBF and gluon fusion results (Eqs. ( p.ip and Eq. ( |2.3D ) we note that 
gluon fusion dominates for small enough or large enough tag jet transverse momenta. At 
the energy scales interesting for LHC measurements there is, however, a very important in- 
termediate k_i_ range where the WBF process dominates. We also note that [M^gl^ depends 
on the azimuthal angle as cos^ whereas \Mzz\'^ and \Mww\'^ are flat in azimuth. 



3. Gluon resummation 



To compute the corrections to the lowest order amplitudes of the previous section arising 
from QCD radiation we follow closely the calculation of the "gaps between jets" process 
presented in [23-28]. This process has been widely studied in the literature [24,25,28,43-47] 
and in experiments at HERA [48,49] and the Tevatron [50,51]. In the gaps between jets 
process, one is interested in final states which contain two hard jets at a scale Q with the 
restriction that there be a region (of size Y in rapidity) between the jets into which no 
further jets with pT > Qo can be emitted. This situation is very similar to the Higgs plus 
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two jets process. The only difference is the presence of a third colour neutral particle (the 
Higgs) in the final state and not surprisingly we can carry over much of the formalism. 

In order to account for the presence of a Higgs, we follow the calculation presented for 
the five parton evolution in [52] in conjunction with the result of [24,25]. In this section, 
we focus upon quark-quark scattering. This is the whole story for the WBF process but 
clearly the gluon fusion process can be generated using either quarks, anti-quarks or gluons 
within the colliding protons. The generalization to processes initiated by anti-quarks or 
gluons is straightforward and we present final results for all of the relevant subprocesses. 
In the qq — qqH case, the final result is 



M = exp 



/ Q 

^ ^as{qT)—T I Mo, (3.1) 



, vrj qr 

\ Qo 



where 

^ Cf 




and 



- 2 /'(^'y3,y4)l+ ( ,_ 1 (3-2) 

^ ^ (^ sinhGy^ + y/2) sinhGy^ + y/2)\ 

2V smhdysl - y/2) smh(|y4| - y/2) y 

In the frame in which the incoming quarks collide head on with equal energies, the rapidities 
of the outgoing quarks jets are 2/3 and 2/4 and Y{< |y3 — y/^\ — 2R) is the region in rapidity 
between the two tag jets into which radiation above Qo is forbidden (i? is the radius 
parameter). The effect of including the Higgs boson is essentially kinematical since, unlike 
the gaps between jets case, we do not require that = — ^4. Numerically the effect is 
small in the eikonal region since /o is a small constant for large enough Y and R not much 
smaller than 1. 

We need also to explain the colour basis implicit in the above. We use a singlet-octet 
basis such that the lowest order amplitude for neutral current quark scattering is 

To reconstruct the cross-section, we must restore the colour factor for the gluon fusion 
process, i.e. we need a = M^SM where 



1 



(3.5) 



The amplitude for charged current fusion, Mww, does not interfere with the neutral current 
fusion processes and soft gluon evolution therefore proceeds as in Eq. with = 

{M\Y]y,0). The amplitudes for qq qqH, qg qgH and gg — > ggH are obtained 
analogously and the relevant evolution matrices are presented in the appendix. 

It is now interesting to reconsider the various radiation models discussed in the intro- 
duction and look at them in the light of Eqs. (|3.1|)-(|3^. Clearly the truncated shower 
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approximation has no direct correspondence in terms of a F. Exponentiating the differ- 
ential probability for emitting an extra gluon, as is done in the soft gluon exponentiation 
model, should catch the diagonal Cp dependent parts of F, as well as the Y dependent 
part in the octet-octet entry. The imaginary parts, coming from Coulomb gluon exchanges 
which generate on shell intermediate quarks, will however not appear when only real emis- 
sions are considered and the octet-singlet mixing parts of F are thereby lost. In Monte 
Carlo parton showers we anticipate that the leading and diagonal Cp dependent terms 
are included. Once again the Coulomb gluon in terms are not included. 

Let us make a few general remarks about the soft gluon evolution. Firstly, the Coulomb 
gluon contributions which generate the in terms are just the exchanges which in QED 
would generate an overall phase in the amplitude whereupon they would have no effect 
on observables. In the non-abelian QCD case they do not, however, simply build a phase 
and they play a potentially crucial role. In particular, it is through these terms that the 
singlet and octet sectors are able to communicate with each other. Indeed Coulomb gluon 
exchange is in this way seen to be responsible for generating the QCD pomeron. 

For large enough Y we may neglect the Cpp dependent terms in F and we refer to this 
approximation as the vr^DLLA. In this approximation, the lowest order in ag where soft 
gluon effects enter into WBF contribution to the cross-section is ©(af). This is not the 
case for octet exchange, where the leading contribution at large Y is O^agNY). Thus we 
expect that the Cpp dependent terms should be relatively more important in WBF than 
in gluon fusion. Note also that the interference between WBF and gluon fusion starts at 
0(0^) in the cross-section. 

Recently Mert Aybat, Dixon and Sterman [36] have confirmed that the two- loop soft 
gluon evolution matrix is proportional to the one-loop matrix F and hence that its effect 
can be accounted for if one uses the strong coupling in the CMW scheme [53], i.e. our 
calculations can be extended to include their result via a re-definition of 



with 




The effect is thus a small downward shift in ^. Although the two-loop evolution matrix is 
an important component of the higher order calculation it is not the whole story and so at 
present a full next-to-leading logarithmic calculation is not possible. 

4. Numerical results 

In Fig. ^ we show the square of the helicity averaged lowest order amplitudes, Eq. ( |2.1| ) and 
Eq. (|2.3| ), in the case of a Higgs produced with zero transverse momentum, i.e. ki_j_ = —k2±- 
We have used a one-loop running with nj = 5, as{mz = 91.2 GeV) = 0.120, a constant 
OEM of 1/127.9 and sin^ Oyy = 0.231. For ZZ fusion we have chosen to display the squared 
amplitude for u-type quarks only, but clearly the result for other partonic combinations 
can be obtained by a simple re-scaling of the displayed curve. 
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In order to obtain realistic hadron 
level cross-sections the parton distribu- 
tion functions must be folded in. How- 
ever we do not perform that step here, 
preferring instead to concentrate on 
the physics associated with soft gluon 
radiation. Saving the interference for 
later, we now move on to consider the 
QCD evolution of the gluon-gluon fu- 
sion and WBF amplitudes. 

The various squared amplitudes, 
divided by the corresponding lowest 
order results, are displayed in Fig. ^ 
as functions of ^ for various values of 
Y and for |y3| = \yi\ = Y/2 + R with a 
cone radius R = 1. To get a handle on 
the full LLA as compared to the vr^DLLA, the full LLA evolution is displayed with solid 
lines and the tt^DLLA with dashed. We can see clearly the effect of the larger emission 
probability for the gluon fusion processes and that the vr-^DLLA typically provides a good 
approximation. Perhaps most striking is the fact that the suppression with increasing ^ 
is much the same for all of the gluon fusion processes, i.e. it is not very much dependent 
upon whether the primary jets are quarks or gluons'^. This is to be expected and is a result 
of colour coherence. Large angle soft gluon emission is not sensitive to the colour charge 
of the external partons, rather it is as if the radiation is from the exchanged t-channel 
gluons. In fact, in the large approximation, the suppression is exactly equal for all 
partonic subprocesses as a result of the diagonal nature of the evolution matrices (listed in 
the appendix), i.e. 

a f 2NY^\ 
— ~ exp . 

O-Q V ^ / 

Collinear gluon emission is a different matter however, and its effect appears through the 
diagonal evolution proportional to p. In that case, radiation off gluon and quark jets occurs 
in proportion to their respective colour charges. For large enough cone radii, the gap is 
sufficiently removed from the direction of the outgoing partons that collinear evolution is 
unimportant (i.e. p is small). Hence the relatively small difference between the solid and 
dashed curves in all of the panes in Fig. |3|. In order to enhance the differences between 
quark and gluon jets, one would need to take a small cone radius in order to enhance the 
role of p ~ log(l/i?) at small R. 

We note that the WBF cross-section (the uppermost curves in the first row of the 
figure) actually rises slightly as we go from y = 3 to y = 6. This is a result of the fact 
that for singlet exchange leading logarithmic soft gluon effects are contained entirely in 
the exp(— ^^/o) term at large Y whilst at lower values of Y the mixing into octet and 
the associated suppression plays a role. Thus at large Y we see the dominance of pure 
singlet evolution in the WBF process. This result is not so transparent in the octet-singlet 

^Of course the radiation in qq and qq processes is very similar. 
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Figure 2: 



The relative squared amplitudes of 



Eq. (2.1) and Eq. (2.3). The solid line is gluon-gluon 
fusion, the long dashed line ZZ fusion and the short 
dashed line WW fusion. 
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Figure 3: The cross-section after soft gluon evolution divided by the corresponding lowest order 
cross-section as a function of ^ for Y ^ 1, Y = 3 and Y ~ 6. The first row corresponds to quark- 
quark scattering, the second to quark-antiquark scattering, the third to gluon-quark scattering and 
the fourth to gluon-gluon scattering. The upper family of four curves in the qq and qq plots are for 
WBF whereas all other curves are for gluon fusion. In all cases, the dashed curves correspond to 
the TT^DLLA and the solid lines correspond to the LLA. Grey curves are the order as expansions 
of the corresponding black curves. 

basis we use here but it is apparent in a basis in which the evolution is diagonal [47]. 
We caution however that the a'^Y"' terms will be important for large Y and (apart from 
the double logarithmic Q;"y" log"'(Q/Qo) subset) they are not included in the approach 
presented here. 

What should also be noted is the importance of soft gluon resummation, especially 
for the gluon fusion process. The region around ^ ~ 0.3 is typical of values which will 
be probed at the LHC, and for such values of ^ the fixed order calculation clearly is 
inadequate. As for the Monte Carlos, the situation is less clear since they do go beyond 
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Figure 4: The cross-section after soft gluon evolution divided by the corresponding lowest order 
cross-section as a function of ^ for F = 3. The solid curves are the full result (as in Fig. ^ and the 
grey curves are the result of ignoring the off-diagonal entries in F. In the qq and qq cases, the upper 
long dashed curve is the colour singlet contribution to WBF, the short dashed curve is the octet 
contribution to the gluon fusion cross-section, the lower long dashed curve is the singlet part of the 
gluon fusion cross-section. For the qg and gg cases, the short dashed line is again the antisymmetric 
octet contribution, the long dashed line is the singlet contribution whilst the medium dashed line 
corresponds to the small symmetric octet contribution. The 10 and 27 exchange contributions 
formally present in the gg case are too small to be visible. 



fixed order. Nevertheless, we know that both the soft gluon exponentiation model and 
the parton showering as manifest in HERWIG and PYTHIA will fail whenever the off- 
diagonal (Coulomb gluon) parts of T become important. The influence of the off-diagonal 
(i.e. octet-singlet mixing) terms is shown in Fig. ^ where the grey curves are obtained by 
neglecting them. Their effect starts to become important for ^ ^ 0.3 (in the case of gluon 
fusion this value reduces slightly as Y increases). 

The colour mixing is related to another striking feature of Fig. ^, namely the fact 
that the evolution of the gluon fusion contribution starts out as an exponential decay, but 
then evolves to a much slower decay. The exponential decay for low ^ is expected since 
both the p and the Y terms in T represent exponential decays in ^. The non-exponential 
property is thus a consequence of the mixing induced by the off-diagonal imaginary parts 
of the evolution matrix. This is clarified in Fig. ^ where the different colour components 
are shown separately. In fact this behaviour is arising because a part of the perturbative 
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QCD pomeron is present in our calculations and it is visible in Fig. |^ as the singlet part 
of the initial octet exchange (long dashed line). For sufficiently large Y the resummation 
we perform here will not be adequate and the BFKL resummation of the series in OgY 
will become important [43,54,55]. The problem of summing simultaneously the soft and 
high-energy logarithms has been explored in [47] where it was illustrated that BFKL effects 
may well be important for Y values as small as 3. In the qq and qq cases, we note that the 
long dashed curve represents not only the singlet contribution to the gluon fusion cross- 
section but also the octet contribution to WBF. That these two are identical in the leading 
logarithmic approximation follows from the symmetry of the evolution matrix in a basis 
where S = diag(l, 1) [56]. 

Finally we turn our attention to the interference between the ZZ fusion and gluon- 
gluon fusion contributions. In our calculations there is no interference if one integrates over 
the cos (j) dependence of the tag jets since soft gluons do not change the original azimuthal 
dependence of the lowest order matrix elements. A sizeable interference would certainly be 
of interest since the azimuthal angular distribution ~ cos (p differs both from the flat WBF 
distribution and the gluon-gluon fusion distribution ~ cos'^ (j). The three different angular 
distributions would in principle provide a handle which would help isolate the separate 
contributions. 

However, we find that the level of interference is very low. In fact, even if the color 
structure is neglected, there is a cancellation between quarks of different helicity and flavour 
already at tree level. In the approximation of vanishing sin^ 6w and assuming only u-type 
quarks, the level of cancellation due to the helicity structure is 1/4, since in only one case 
out of four will there be two left handed quarks interacting. This reduces to 1/36 if we 
assume that the proton is made up of two up quarks and one down quark which in turn 
reduces dramatically to 1.6 x 10"'^ when we restore the correct value of the weak mixing 
angle. ^ 

Writing the interference between ZZ and gg fusion contributions as 

/ = M^SM-M+^SMgg -M^^SM^z, (4.1) 

where M = M^^ -|- M^^ is given in Eq. (^T]), we show —I/gq in Fig. ^ for the case of u 
quark scattering: as anticipated it is very small. 

5. Conclusion 

We have performed a partial summation of the leading soft gluon contributions in Higgs 
plus two jet production in the presence of a (mini-)jet veto. The imposition of such a 
veto helps enhance the contribution from weak boson fusion by exploiting the fact that 
no colour is exchanged between the outgoing tag jets (and hence radiation is suppressed 
in the central region). We flnd that large differences in (mini-)jet activity can indeed be 
expected depending on whether the hard process is a color octet gluon fusion process or a 
color singlet weak boson fusion process. 

^In fact the interference would have vanished identically if sin^ 9w = 1 /4. 
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Figure 5: The size of tlie ggH-ZZH interference (normalized to the total lowest order cross- 
section for gg and ZZ fusion) in the case of u-type quarks only. In (a) F = 1, (b) F = 3, (c) 
Y = 6. Solid curves are for |fc_L| — 30, long dashed curves for |fej_| = 100 and short dashed curves 
for Ifc, I = 300. 



Neglecting non-eikonal contributions (which is a good approximation when the tag jets 
are at low angles relative to the collision axis) the probability for extra gluon emission in 
the gluon fusion case is proportional to agNY and hence increases as the gap between the 
tag jets increases. For the electroweak case it is proportional to agCpp which implies that 
radiation is not strongly enhanced as the gap increases. This is apparent in the vr^DLLA 
where one keeps only the leading Y and ivr dependent terms. In that case, corrections to 
the WBF process starts at ©(af). Apart from containing the expected general features, 
our results also indicate when it is necessary to go beyond fixed order perturbation theory 
and resum the logarithms in the ratio of the tag scale to the veto scale. We conclude that 
for both the gluon fusion and weak boson fusion production of a Higgs, resummation is 
important in the region which can be probed at the LHC. Moreover, our results indicate 
that for large enough values of the evolution variable, ^, it is not sufficient to account for 
the resummation using existing parton shower algorithms or the soft gluon exponentiation 
model. 

We have also made an estimate of the size of the interference contribution between 
weak boson fusion and gluon-gluon fusion. The result is very small, never exceeding 1% of 
the lowest order cross-section. 
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A. Evolution matrices 

In this appendix we list the evolution matrices for the various subprocesses relevant to 
Higgs production through gluon fusion. The function p{Y,ys, 1/4) is as defined in Eq.(^^) 
and this is the only difference compared to the evolution matrices presented in [23,25,28]. 
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For qq qqH the evolution proceeds in a similar fashion to the qq case except that 
the evolution matrix, in the singlet-octet basis, is 



Cjp I 

^^(^'^3, 2/4)1+ . 
2 \ -ZTT ( 



2 2JV 



^7^ 



(A.l) 



For — > gg'ii" we use the basis 



ci 



C2 = C^rsr2c(r'')^^^^ 
C3 = ifrBr2c{T'')^^^^ 



(A.2) 



corresponding to singlet, symmetric octet and antisymmetric octet t-channel exchange. In 



this case Sg^ is 



/ N{N^ - 1) 
S..= ^(iv2-4)(iV2-l 

\ 

and the evolution matrix is 

/ N-Ki 









J 



■ 19 



(3Ar2-i) 1 







2TTi \ 
2NY Nni 



(N'^-i) 

\ Am ^ ^ > m 2NY ) 



For gg — ggH we use the basis (setting A'^ = 3) 
3 



C'a — ^frAricfrBr2Ci 

C4 = — (<JrArs'^rir2 



2 (^'■A»'s"''ir2 



(^n ro ~l~ /t ro (^i 



gOj-^r-^O^grj ^"'rAric"'rBr2C- 



In which case 



"99 



/l \ 
8 
8 
20 

yo 27/ 



and 



^99 = 2p{Y,yz,yA) 1 + 



/-fvri 

2^ 

3 







-37ri 








-Im 3r-f7ri -^7ri 
V -^TTi -|7ri 4y-|7rz/ 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 
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